Quantum phase slip interference device based on superconducting nanowire 
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We propose a transistor-like circuit including two serially connected segments of a narrow super- 
conducting nanowire joint by a wider segment with a capacitively coupled gate in between. This 
circuit is made of amorphous NbSi film and embedded in a network of on-chip Cr microresistors en- 
suring a sufficiently high external electromagnetic impedance. Assuming a virtual regime of quantum 
phase slips (QPS) in two narrow segments of the wire, leading to quantum interference of voltages 
on these segments, this circuit is dual to the dc SQUID. Our samples demonstrated appreciable 
Coulomb blockade voltage (analog of critical current of the SQUIDs) and periodic modulation of 
this blockade by an electrostatic gate (analog of flux modulation in the SQUIDs). The model of this 
QPS transistor is discussed. 
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The discovery of the Josephson effect [ij triggered the 
development of superconductive electronic devices oper- 
ating on the principle of the classical behavior of the 
collective quantum variable, i.e. the superconducting 
phase difference. These include rf and dc SQUIDs 
single-flux-quantum logic circuits Josephson voltage 
standard arrays 4], microwave receivers [s!], etc. The op- 
eration of these devices is based on the nonlinear {2-k 
periodic) dependence of the supercurrent I = Ic sin (p 
on the phase difference ^ on the Josephson junction 
(JJ), where Ic is the critical current. The electromag- 
netic impedance of the JJ for a small harmonic sig- 
nal is therefore expressed in terms of Josephson induc- 
tance Lj which varies with ip also in a periodic fashion, 
L^^{ip) = (27r/c/<I>o) cos where <I>o = h/2e is the flux 
quantum. 

Recently, a new type of nonlinear superconducting ele- 
ment, dual to the JJ, - the quantum phase slip (QPS) 
nanowire (see, for example, the review on the phj^s- 
ical properties of superconducting nanowires Ref. 
and references therein) - was proposed in Refs. 0, 
Being included in a superconducting ring with flux- 
bias near degeneracy points, $(, = (m + 0.5)$oi 
is an integer number, the QPS element enables cou- 
pling of quantum states |m) and |m -I- 1). This leads 
to anticrossing of corresponding levels with the gap 
AE = Eqps- For a uniform nanowire of length I, 
the characteristic coupling energy i^qps = ^es with 
es ~ [kBTc/OiRQ/R'O^^^/eMO-SaRQ/R'O, where 
is the critical temperature, ^ the coherence length, R' the 
normal-state resistance per unit length, Rq = hjAe^ ~ 
6.5 kil, the resistance quantum, and a ~ 1 @,J3]- The 
two-level spectrum of this circuit was proposed [7!] for us- 
ing in a flux qubit (compare with the dual charge qubit 
in which the coupling of charge states is possible due to 
finite Josephson energy Ej — {<^q/2ti)Ic [^)- 

Coupling of the nanowire to a high-ohmic environ- 
ment through the appropriate boundary conditions im- 
poses dramatic changes in its dynamics [l3]- Based 
on the quantum duality of the Josephson and QPS ef- 



fects stemming from conjugation of the phase and charge 
variables, Mooij and Nazarov had further predicted [8| 
that in the charge-bias regime, the QPS element should 
demonstrate the Coulomb blockade behavior with pe- 
riodic dependence of voltage on the injected charge 
viz. VQPs(g) = V"c sin(7rg/e) with Vc = 27r£;QPs/2e. 
Thus, this element should behave as a nonlinear capac- 
itor, Cqpg((7) = rfVQps/rf? = {ttVc/ e) cos{-Kq/ e) with a 
2e-periodic dependence on q. 

It is clear that such a QPS element could be exploited 
in the engineering of electronic devices controlled by the 
charge in the classical manner. These devices are dual 
to the Josephson ones controlled by the classically be- 
haved phase. For example, applying microwave irradia- 
tion of frequency / should lead to the formation of the 
current steps in the I-V curve at / = 2nef (n is in- 
teger), enabling a fundamental standard of current [1]. 
These steps are dual to Shapiro voltage steps in a JJ I-V 
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FIG. 1: (Golor online) (a) The layout and (b) the simpli- 
fied electric circuit diagram of the QPS transistor embedded 
in the 4-terminal network of on-chip resistors. The device 
includes two QPS elements denoted by diamond symbols, ki- 
netic inductances of the nanowire segments and a capacitive 
gate. Thicknesses of NbSi and Gr films are 10 nm and 30 nm, 
respectively. 
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curve a,t V = n$o/ 111- Developing such an analogy, in 
this paper we proposed and realized a QPS-based single 
charge transistor which is dual to dc SQUID and can be 
operated as an electrometer. 

Our device (see schematic diagram in Fig.[T^) has a 
symmetric in-line configuration and includes two nar- 
row pieces of superconducting nanowire joined by an is- 
land (a wider middle part of the same nanowire) with 
a capacitively coupled gate. This transistor is embed- 
ded in the network of compact high-ohmic resistors, 
R ^ 0.4 Mfi 3> Rq, which ensure a sufficiently high elec- 
tromagnetic impedance seen by the transistor. Due to 
this improvement, a classical charge regime of operation 
with significant damping of dynamics is realized. (Note 
that the recent proposal by Hriscu and Nazarov [13] deals 
with a complementary QPS transistor, without resistors, 
operating in a quantum regime, which, in contrast to our 
circuit is phase biased and, therefore, not strictly dual to 
that of the dc SQUID.) 

The Kirchhoff equation for the equivalent electric cir- 
cuit presented in Fig.[T)D has the form 

V1+V2 = Vb, = Lk^q^+R^q^ + VQPS^{q^), i = l,2. (1) 

Here, VQpsi(<Zi) = VciSm{TTqi/e) are the 2e-periodic volt- 
ages on the QPS elements on either side of the island. 
The charge conservation relation takes the form 



qi - 92 



iCg + Co){Vi-V2). 



(2) 



This equation refiects the balance of the charges injected 
in the island and the polarization charges induced by the 
gate and fioating potential of the island having a self- 
capacitance Co with respect to ground. In the case of 
sufficiently small gate- and self-capacitances. 



(3) 



and small injected currents (viz. ±51^2), the difference 
charge gi — 52 — —CgVg = Qg is totally controlled by the 
gate voltage Vg. 

One can easily find that equations Eqs. ([TJS]) are dual 
to those of the dc SQUID, see, for example, Ref. [13]. The 
charges qi and q2 are dual to the phases across individ- 
ual JJs, whereas amplitudes Vci,c2 play the roles of the 
critical currents. To continue this analogy we can put in 
correspondence the pairs of parameters Lki^k2 and JJ ca- 
pacitances, Ri^2 and shunting resistors, etc., for the QPS 
transistor and dc SQUID, respectively. The condition of 
small gate- and self-capacitances Eq. ([3]) is dual to the 
requirement of a small inductance of the SQUID loop 
(L ^ Lj). Similar to external flux control of quantum 
interference in the SQUID pji], the voltage on the QPS 
transistor is resulted from quantum interference of volt- 
ages on individual QPS elements and controlled by the 
gate charge Qg, 



where Q = qi + erj/ir = q2 + erj/n — Qg is the "average" 
charge. 



tanry 



2Vc2 tan(7rgg/2e) 



(5) 



Vc+ + Vc- tan2(7rQg/2e) ' 

Vc± = Vci ± Vc2 , and the Coulomb blockade voltage 

= Vol + Vc2 + 2K1K2 cos{nQg/e). (6) 

In the case of a symmetrical circuit, Vd = Vc2 = Vc, Q ~ 
qi 2 ± Qg/^, Eq. ^ yields the maximum modulation of 
the blockade voltage, Vm = 214] cos(7rQg/2e)]. In the case 
of a highly asymmetrical circuit, say Vc2/Vci — a 1, 
charge Q k, qi and the blockade voltage is only slightly 
modulated by the gate charge Qg, 



Vm = Vci[l + acos{TrQg/e)]. 



(7) 



Thus the equation of motion Eq. ([T|) takes a form of the 
resistively shunted junction (RSJ) model for dc SQUID, 



LkQ + RQ + V,n{Qg) sin(^Q/e) = Vb 



(8) 



with Lk = Lki + Lk2 and R = R2 + R2- 

Equation (|5]) describes the dynamics of a nonlinear 
oscillator with finite damping. It yields a dc jy-curve 
with both the static- and running-charge (the oscillating 
regime, uj = Tr{Q)/e = 7r(/)/e) branches, dependent on 
Qg in a periodic fashion. The dimensionless parameter, 
crucial for the dynamics, introduced earlier in Ref. Q, 



^cLk/R, 



(9) 



where Wc = ■nVm{Qg)/eR is a characteristic circular fre- 
quency. Eventually, /3qps is the analog of the Stewart- 
McCumber parameter in the Josephson dynamics. In 
the most realistic case of large damping, /3qps < 1, the 
/F-curve has a shape with characteristic back bending. 



{v)^[v^-SQ,) + R^{i?Y'^-m, 



(10) 



V = Vqpsi + Vqps2 = VyniQg) sin(7rQ/e) 



(4) 



where (y) is the average voltage on the nanowire. 

The fabrication method (some modification of the 
method developed in Ref. jlj]) included the following 
steps. The samples were fabricated on a silicon substrate, 
having 300 nm thick thermal oxide layer, combining two 
processes: shadow evaporation and the sputter deposi- 
tion/etching step. The Cr resistors were fabricated in a 
single vacuum cycle with the AuPd contact wires and 
micropads utilizing the shadow evaporation technique 
through a bi-layer PMMA/Copolymer stencil mask. The 
30 nm thick Cr resistors were evaporated first at a low 
residual pressure of oxygen (^ 10~^mbar) followed by a 
50 nm thick layer of AuPd from an angle for which the 
narrow stencil openings for Cr resistors were overshad- 
owed by the mask. Using this trick, the formation of 
AuPd shadows parallel to the Cr resistors was avoided. 



3 




-2-101234 

Voltage (V) (mV) 



FIG. 2: (Color online) The charge-modulated IV-curves of 
Sample A recorded in a current bias regime for two gate volt- 
ages shifted by a half-period. In the region of small currents 
(blown up in inset) one can see the modulation with period 
AJ = 13.5 pA, which is due to the asymmetry of off-chip bias- 
ing circuitry resulting in the current dependence of the electric 
potential of the transistor island, 6V = (-Rbiasi —Rbia.s2)I and, 
therefore, of the effective gate charge, 5Qg = (Cg -I- Co) 5V (cf 
Eq. (d}). The green dashed line shows the shape of the bare 
IV -crave given by the RSJ model Eq. (|10[) . 
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FIG. 3: (Color online) The /T/-curves of Sample B measured 
in the voltage bias regime at different values of gate voltage 
Vg. The bottom right inset shows details of the Coulomb 
blockade corner. Upper left inset: the gate voltage depen- 
dence of the transistor current measured at different bias volt- 
ages Vf, providing a steady increase of (V) from 0.321 mV up 
to 0.481 mV in 4 mV steps (from bottom to top). 



The Cr resistors and other parts of the circuit were pro- 
tected by a PMMA mask while 10 nm thick amorphous 
NbajSii-a; film (x — 0.45, defined by calibrating the sput- 
ter rate for each element at a given sputter power and 
periodically confirmed by EDX measurement using sepa- 
rately deposited films on Ge substrates) was co-sputtered 
on the substrate making contact only with the AuPd mi- 
cropads. These micropads allowed making an electrical 
contact of the NbSi film with the Cr resistor circuitry. 
To remove organic residue and water for the purpose of 
forming reliable contacts, we cleaned the micropads in 
RIE with oxygen plasma and baked the sample overnight 
in N2 atmosphere at a temperature of 120 C° before NbSi 
deposition. The subtrate was rotated 20 times during a 
one-minute-long deposition. After the lift-off, the wafer 
was coated with inorganic negative tone HSQ resist (XR- 
1541, Dow Corning) patterned with e-beam. After expo- 
sure, the HSQ resist had characteristics of thermal silicon 
oxide and could be used as an etch mask to define the 
NbSi nanowires along with the island and the gate struc- 
ture made out of the same NbSi film. In a final step, the 
ICP etch process with SFg gas was used for etching. 

The samples were measured in a dilution refrigera- 
tor with electrical lines equipped with microwave fre- 
quency filters made of pieces of Thermocoax"""^ cable. 
We used battery power sources and home-made electron- 
ics for either a voltage or a current bias of our circuits. 
Most of the measurements were performed at the low- 



est temperature of the fridge, T = 15 mK. The layout 
allowed independent characterization of both pairs of 
on-chip resistors. Their /V-curves at millikelvin tem- 
peratures were practically linear with resistance about 
20% above the values measured at room temperature. 
The NbSi films had the superconducting transition tem- 
perature Tc « 1 K and the normal state resistance per 
square about 550 il (approaching the level of ca. 600 17 
of the thickness-driven superconductor-to- insulator tran- 
sition in somewhat thicker films of Nb^jSii-^r with sto- 
ichiometry x = 0.15 [isj). yielding for the narrow seg- 
ments the specific resistance R' « Slkil/ fim. The NbSi 
film parameters were measured at T = 3.4 K on the 
stand-alone test wires 20 nm in width, fabricated on the 
same chips. 

Here we report on two samples (A and B, see their char- 
acteristics in Figs. 2 and 3, respectively). In comparison 
with Sample B whose dimensions are shown in Fig. la. 
Sample A had a twice longer island (25 /im) and some- 
what shorter QPS sections (4.5 /.tm each). The layout of 
the gate was identical in both samples. Both samples 
showed an appreciable Coulomb blockade, (/) = 0, and 
a gate effect at small bias. An increase in current (/) > 
1 nA was accompanied by a steep rise in voltage {V) indi- 
cating a gradual turning of the wire into the normal state. 
On a large scale, the /F-curves of both samples exhibited 
a positive excess current, i.e. had a shifted linear asymp- 
tote, /off - (/) - {V)id{I)/d{V))i^v)^^ « 1-5 - 2.0 nA, 
indicating that part of the current at large bias was still 
transported by a superconducting component [16]. This 
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shape of the IV-cmve argues that the observed gate ef- 
fect has a superconductive origin and ahows us to rule 
out the strikingly different effect of the Coulomb block- 
ade in a normal single electron transistor, characterized 
by a pronounced voltage offset [17| or, equivalent ly, a 
negative "excess" current. 

On a small scale, one can see the clear gate mod- 
ulation of IV -curves, although with very different pe- 
riods, AVg ~ 2mV and 150 mV, for Samples A and 
B, respectively. Assuming the 2e-periodic charge de- 
pendence in Sample A (see inset to Fig. 2) capacitance 



Co 



2e/AVg takes the value of about 160 aF. This 



value is not far from the modelled value of 200 aF. The 
depth of the modulation of about 40% suggests - accord- 
ing to our model Eq. ([7]) - the ratio of QPS energies, 
a = Vc2/Vci = £'qps2/£^qpsi ~ 0.2 whereas the average 
value is about 20 /ieV. Inserting this value in the formula 
for the QPS energy and assuming the length of the QPS 
junction of the order of the coherence length ^, we obtain 
the ratio Rq/R'£, w 9 and, therefore, a rough estimate 
^ « 20 nm. This value is close to the values ^ = 10—15 nm 
found in Ref. Q. 

The shape of a bare (not modulated) IV-cmve is only 
qualitatively similar to that given by Eq. (jlOp and shown 
by the dashed line. A back bending weaker than in theory 
can be attributed to the effect of stray capacitance of the 
resistors, resulting in a roll-off of effective impedance R 
in Eq. (|10l) with the rise of current (oc ui). Rounding of 
the Coulomb blockade corner is attributed to the effect 
of noise, omitted in our model. 

Sample B demonstrated an even stronger gate effect 
(see upper left inset of Fig.[3l), but also the peculiar 
properties. A huge period AVg « 150 mV of the gate 
voltage dependence corresponds to coupling capacitance 
Cg « 2 aF. This value is two orders smaller than the de- 
signed value, so we conclude that the actual island in 
this sample had much smaller dimensions (very crudely, 
on the order of 100 nm) and located presumably in one 
of the narrower segments of the nanowire. The available 
single gate did not allow any conclusions about exact size 
and location of this island, whereas the depth of modu- 
lation of about 20% yielded the ratio of QPS energies, 
-Bqps2/-Eqpsi « 0.1 with E'qpsi ~ 150 neV. A smah 
hysteresis and almost vanishing back bent part in the 
IV-cmve, seen in the lower bottom inset of Fig.|3l may 
indicate that the value of the damping parameter /3qps 
Eq. ([9|) is between 1 and 2 [l3| , whereas our estimation 
yields a value of /3qps smaller by an order of magnitude, 
ensuring a heavily overdamped regime. 

Interestingly, the observed periodic pattern in Sam- 
ple B was superimposed on another one having a sub- 
stantially larger period and a much weaker modulation, 
discernible in the traces at smaller bias in the upper 
left inset of Fig.[3l This behavior can be due to the 
emergence of an additional (very small) island neighbor- 
ing in line with the actual island. Thus, in spite of a 



rather high homogeneity in thickness and width, the en- 
tire nanowire can be considered as a circuit including sev- 
eral weaker sections with a local increase in QPS energies 
per unit length, es(a;). In the spirit of the single charge 
interferometer operation, Eqs. (HHHl), the contributions of 
these weak sections to the total QPS energy should be 
summed up, taking into account phases proportional to 
the charges induced by the gate on all intermediate is- 
lands. The resulting quasiperiodic dependencies can be 
also interpreted in terms of the Aharonov-Casher effect 
demonstrated in the experiments with phase-biased ar- 
rays of small JJs (i.e. in a lumped-element analog of 
superconducting nanowire) [lil . Recent calculations 
by Vanevic and Nazarov [20| support the hypothesis that 
apparently homogeneous nanowires may naturally have a 
strongly inhomogeneous distribution of specific QPS en- 
ergy es{x) because of the exponential dependence on the 
local parameters. Large spatial fluctuations of the local 
energy gap detected by scanning tunneling methods were 
also reported for thin disordered films of TiN [21], NbN 
^ and InO 

In conclusion, we have demonstrated the single-charge 
effect in superconducting nanowires having a transistor 
configuration with a capacitively coupled gate, embedded 
in a high-impedance environment. A deeper understand- 
ing and better control of the nanowire parameters deter- 
mining the characteristics of these transistors and other 
possible circuits is urgently needed and motivates us to 
conduct further research. Generally, the demonstrated 
duality of the QPS transistor and the dc SQUID may 
open the way towards interesting applications of QPS 
nanocircuits in electronics and metrology. 
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